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Abstract 

We study the impact of contagion in a network of firms facing credit risk. 
We describe an intensity based model where the homogeneity assumption is 
broken by introducing a random environment that makes it possible to take 
into account the idiosyncratic characteristics of the firms. We shall see that our 
model goes behind the identification of groups of firms that can be considered 
basically exchangeable. Despite this heterogeneity assumption our model has 
the advantage of being totally tractable. The aim is to quantify the losses that a 
bank may suffer in a large credit portfolio. Relying on a large deviation principle 
on the trajectory space of the process, we state a suitable law of large number 
and a central limit theorem useful to study large portfolio losses. Simulation 
results are provided as well as applications to portfolio loss distribution analysis. 
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Introduction 



During the last years the challenging issue of describing the dynamics of the loss pro- 
cess connected with portfolios of many obligors has received more and more attention. 
Applications can be found both for management purposes (see [8j) and in the liter- 
ature dealing with pricing and hedging of risky derivatives as CDOs (Collateralized 
Debt Obligations). For a discussion of this framework see [15] and [18]. 

When dealing with portfolio losses, it becomes crucial the modeling of the depen- 
dence structure among the obligors. One standard procedure is to directly specify the 
intensity of default of the single obligors belonging to the portfolio in order to infer 
the dynamics of the global system and thus the distribution of the aggregate losses. 
In the context of reduced form models a rather general framework is the conditionally 
Markov modeling approach by Frey and Backhaus (see [12] and [I3]). One drawback 
of the intensity based models is the difficulty in managing large heterogeneous port- 
folios because of the presence of many obligors with different specifications. In this 
case it is common practice to assume homogeneity assumptions in order to reduce the 
complexity of the problem. A typical approach is to divide the portfolio into groups 
where the obligors may be considered exchangeable. 

In this paper we describe an intensity based model where the homogeneity as- 
sumption is broken by introducing a random environment that makes it possible to 
take into account the idiosyncratic characteristics of the firms. We shall see that our 
model goes behind the identification of groups of firms that can be considered basically 
exchangeable. Despite this heterogeneity assumption our model has the advantage of 
being totally tractable. 

The goal is to describe the evolution of the losses for a large portfolio where het- 
erogeneity and direct contagion among the firms are taken into account. We denote by 
L^{t) the random variable describing the losses at time t G [0, T] for a portfolio of size 
A^. Our approach works as follows. First we study the N ^ oo limiting distributions 
on the path space of some aggregate variables useful to characterize the evolution of 
L^{t) for t G [0,T]. To this effect we shall derive an appropriate Law of Large Num- 
bers based on a Large Deviations Principle in order to describe a limiting behavior 
that can be considered as a asymptotic regime with infinite firms. Finally, we study 
the finite volume approximations (for finite but large A^) of the limiting distribution 
via a suitable version of the Central Limit Theorem that describes the fiuctuations 
around this limit. In most cases, these dynamical fiuctuation theorems are proved by 
method of weak convergence of processes; this approach has been widely applied to 
models close in spirit to this work. We quote [12j and ^ for applications to finance. 
The effectiveness of those methods for heterogeneous models is, however, unclear. We 
follow here a different approach, which allows to prove a Central Limit Theorem di- 
rectly in the underlying trajectory space. This approach is based on a general Central 
Limit Theorem in P]. Although various applications of this theorem to fluctuations of 
Markov processes can be found in the literature (see e.g. [T] and [S]), to our knowledge 
the first application to a non-reversible Markov process. 



2 



In the risk management context, our model may be useful for the management 
of large portfolios, in the spirit of other models proposed in [16] or in It has 
been remarked that in many real world applications default are rather rare events, 
so that, for instance, the fraction of defaulted firms is close to zero and a normal 
approximations is not meaningful. Our models and results are only concerned with 
time scales for which a proper fraction of the portfolio is likely to be affected by the 
defaults. 

We believe that our paper may be considered as an original contribution in the 
modeling of portfolio losses dynamics that accounts for both heterogeneity and con- 
tagion. On the other hand, to our knowledge, this is the first attempt to apply large 
deviations and normal fluctuation theory on path spaces (that is, in a dynamic fash- 
ion) for finance or credit management purposes, except for what contained in[7]. For 
a survey on large deviations methods applied to finance and credit risk see |17] . 

The models we propose in this paper are the simplest heterogeneous models describ- 
ing systems comprised by many defaultable components, whose defaults are positively 
correlated via an interaction of mean-field type, i.e. with no geometric structure. Al- 
though we have been inspired by financial applications, we believe the basic principles 
should apply to other contexts. 

The outline of this paper is as follows. In Section 1 we illustrate the model and 
the main theorems. In Section 2 we apply these results to the large portfolio losses 
analysis. Some examples with explicit computations and simulations are also provided. 
In Section 3 we conclude. Appendix A is devoted to the proofs of the three main 
theorems stated in Section 1. 



1 Model and main results 

Consider a network of defaultable firms, whose states are denoted by yi, 1/2, ■ ■ ■ , Vn, 
Hi G {0, 1}. The event {yi = 1} means that the i-th firm has defaulted. The values 
Hi, 1/2, ■■■ ,yN give rise to an aggregate variable which indicates the global state of 
the network: 



rriN ■■ = 




1=1 



where ai, 02, • • • , Q^Af are given nonnegative numbers, ai can be interpreted as the 
impact the default of the i-th firm has on the aggregate variable mi\f. In order to 
model contagion, we assume the instantaneous rate of default of the i-th firm is an 
increasing function of mAr. More specifically, we assume the rate of default of the i-th 
firm is given by 

where is the indicator function of the set A, and A > 0, 7^ G M are given constant. 
Pi represents the sensibility of the z-th firm to variations of the aggregate variable mjsf, 
while 7j can be interpreted as the "robustness" of the i-th firm: a large value of 7^ 
means that the i-th firm is very unlikely to default within a given time. 
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Thus, for any fixed value of a := (ai,a2, • • • ,aAr), P '■= {Pi, P2, ■ ■ ■ , Pn) and 7 : = 
(71, 72, . . . , 7Ar), the variable y := (yi, 1/2, • • • , Vn) evolves as a Markov chain in contin- 
uous time, with infinitesimal generator given by 

N 

^aAlfiy) ■■=Y.^{y.=^}e'''^''-'V{t) " f{y)l (1-1) 

i=l 

where denotes the configuration obtained from y by changing yi from to 1. We 
assume the system start at time t = from the configuration ?/(0) = (0, 0, . . . , 0). The 
evolution randomly drives the network towards the trap state (1, 1 . . . , 1), which is 
reached in finite time. 

Prom now on we denote by Aj := (a^. Pi, 7^) the triple of the parameters correspond- 
ing to the z-th firm. The Aj's model the heterogeneity of the system. We consider 
here the point of view of disordered models, i.e. we assume Ai, A2, . . . , Aat to be i.i.d. 
random variables, with a given law /i. In order to avoid inessential difficulties, the law 
IX is assumed to have compact support in IR+ x M+ x M. Note that, for a given the 
random variables ai,Pi,'yi are not assumed to be independent. Sometimes, the vector 
A := (Ai, A2, . . . , Aat) will be referred to as random environment. 

Consider a time T > 0, and denote by y[0,T] = (y(t))te[o,T] the trajectory described 
by the configuration under the stochastic evolution (11. ip . Each component ?/j[0,T] is 
either identically or it fiips from to 1 at the default time 

r, := inf{t > : y,(t) = 1}. (1.2) 

By convention, we set ?/i(rj) = 1. This set of {0, l}-valued trajectories is denoted by 
D[0, T]. Each trajectory in T)[0, T] can be identified with its default time (which is set 
to be equal to T if there is no default); thus D[0,T] inherits the topology induced by 
the usual topology on M for the default time. Equivalently, the topology on X'[0,T] is 
the one induced by the Skorohod topology on the set of M-valued functions which are 
right-continuous and admit limit from the left at any point of [0,T] (see e.g. |llj). 

In this paper we are interested in the asymptotic behavior, as ^ +00, of empirical 
averages of the form 

i=l 

where / : 'D[0, T] ^ M is a Borel measurable function, and 

1 ^ 

i=l 

is called empirical measure. More generally, we shall consider the empirical measure 

1 ^ 

PN{y[0,T],\) := j;^'^Sy^io,TlXo 

i=l 
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which is a random measure on D[0,r] x M^. Note that pN{y[0,T]) is the marginal of 
Piv(#,T],A)onD[0,r]. 

In what follows, we denote by Mi the set of probability measures on !D[0, T] x Supp{fi), 
while M will denote the set of signed measure on !D[0, T] x M^. Both sets are provided 
with the weak topology. 
For y[0,T] G 2)[0,T] with y{T) = 1, we set 

T{y[0,T]):=M{t>0:y{t) = l}. 

For Q G M we define 
F{Q) := f Q{dy[O.T],dX)l ^(1 - (l - e-V/«('^*^]''^^>'''W) cit 



+ yiT) 



1.2) 



-7+ /5 / Qidr][0,T],dX')a'r]{t- 

\ J J t=T{y[0,T])_ 

where A = 7) and A' = 7'). For a fixed A, the infinitesimal generator 

(11. ip . together with the initial condition y(0) = (0, 0, . . . , 0) induces a probability 

on I)^[0,T]. We think of as the conditional law of the process given the random 
environment. We denote by 

PN{dy[0, T], d\) := P^{dy[0, T]) ® /i^^(rfA) 

the joint law of the process and the environment. The distribution of PAr(y[0, T], A) 
under Pn will be denoted by Pat o . 

A special case is when all components of A are zero. In this case each firm defaults 
with rate 1, independently of the others. We denote by W the law on D[0,T] of this 
process. 

In what follows, for Qi, Q2 ^ 3Vti, we denote by 

i^(Qi|Q2):=| ^^^^(S^1°SS) ifQi«Q2andglogg^GLHg.) 
[ +00 otherwise 

the relative entropy of Qi with respect to Q2- 

Theorem 1 . The sequence Pn o p'^ of elements of Mi satisfies a Large Deviation 
Principle (LDP) with good rate function 

I{Q) ■=H{Q\W ®p)-F{Q). 

The proof of Theorem [H as well as of the other results stated in this section, is 
postponed to the appendix. 

We recall that the above statement means that, for each Borel subset A of Mi, 
- inf I{Q) < liminf ^ logP^v o p]^\A) < limsup ^ logP^v o pJf^A) < - mi_I{Q), 

o 

where A and A denote the interior and the closure of A respectively; moreover the 
function /(■) is nonnegative, lower-semicontinuous, and the level sets {Q : I{Q) < 1} 
are compact, for each / > 0. 
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Theorem 2 . The equation I{Q) = has a unique solution Q^, that can be identified 
as follows. Consider the nonlinear integro- differential equation 

r |gi(A) = e-''exp[(3jf,{dX')a'qt{X')]{l-qt{X)) 

I go(A) = ^^■^> 

for a real-valued qt{X), t > 0, X = (a,/?, 7) G M.^. This equation has a unique solution 
< qt{X) < 1. For every X fixed, consider the Markov chain on {0,1} with time- 
dependent infinitesimal generator 

Lff{s):=cf{sMl-s)-f{s)], (1.4) 

where 

cf{s) := (1 — s)e~'^exp 

and starting from s = (note that this process jumps only once, from to 1, and it is 
then trapped in 1). Let he the law of this process on !D[0,T]. Then 

Q* = Ql® /i- 

Moreover 

qtiX) = Ql{y{t) = 1). 

Theorems [1] and [2] have a simple consequence. Let U be an open neighborhood of Q* 
in Ml. By Theorem [2], lower semicontinuity of /(■) and compactness of its level sets, 
a standard argument shows that k{U) := miQ^jj I{Q) > 0. By the upper bound in 
Theorem [T] there exists C > such that 

PNiPN eU)< Ce-^^(^), 

thus converges to zero with exponential rate. We summarize this fact in the following 
law of large numbers. 

Corollary 1 . Let d{-, ■) be any metric that induces the weak topology on Mi. Then 
for every e > 0, the probability 

PN{d{pN,Q*) > e) 
converges to zero with exponential rate in N. 

Next result is about the fluctuations of pat about Q*, which has the form of a Central 
Limit Theorem. In most cases, these dynamical fluctuation theorems are proved by 
method of weak convergence of processes. A typical tool in this context is Theorem 
1.6.1 in [llj; it has been widely applied to models close in spirit to this work, see for 
instance ^ and [7j. The effectiveness of those methods for heterogeneous models is, 
however, unclear. The main point is that via Theorem 1.6.1 in [XT] one obtains the 
dynamics of the fluctuation process, which is infinite dimensional; to get a computable 



(3 / fi{dX')a'qt{X') 
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expression for the asymptotic variance of a given function of the trajectory may be 
not feasible. 

We follow here a different approach, which allows to prove a Central Limit Theorem 
directly in the space M. This is inspired by the seminal work of E. Bolthausen |3] , and 
it has been carried out in a context similar to ours in pLj and The main difference 
here is that the underlying stochastic dynamics are not reversible; the related 
difficulties have forced us to introduce a further assumption, that we call reciprocity 
condition: 

(R) There exists a deterministic 6 > such that for alH > 1 the identity (3i = hai 
holds almost surely. 

In economical terms, this means that the sensibility of a firm to variation of the 
aggregate variable m^r is proportional to the impact the default of that firm has in the 
network. In different words, the interaction between firms is symmetric: if the i-th 
firm strongly interact with the network (i.e. is large) then it has a large influence 
on the network but, symmetrically, it is also strongly influenced by the state of the 
network. This assumption may be reasonable in many situations. 
From now on, whenever condition (R) is assumed, (Aj) will denote the pair (aj,7j); 
accordingly. Mi and M will be spaces of measures on !D[0,T] x Supp{p), where 
Suppi/j.) C M+ X M. 

In order to state our Central Limit Theorem, we need to introduce some notations. 
Let Mo be the subset of M comprised by the signed measures with zero total mass. 
Let be the law, induced by Q*, of the Mq- valued random variable 5(y[o,T],A) — Q*- 
We then denote by C;, the space of bounded, continuous, real valued functions on 
D[0, T] X R^. For G Cfe, define G Mq by 

4>{A) := j u,{dR) (^R{A) j (PdRj , (1.5) 

for A C D[0,T] X M2 measurable. 

Theorem 3 . Let (pi, 02, . . . , 0n £ Cfe- Then the Pn-Iclw of the random vector 

^ ( / (pidpN - I (pidQ^ 

J i=\ 

converges weakly as N +oo to an n-dimensional Gaussian probability measure with 
zero mean and covariance matrix C = (Cjj)"^^]^ given by 



Ci,, := j (0. - 0*)(0j- - <l>*j)dQ* - D'F{Q,)[<P,, 0,-], 

where 0* := j (pidQ^, and £)^F((5*) [0j, 0j] is the second Frechet directional derivative 
of F at in the directions 



DF{Q, + h4>j)[4>4~DF{Q,)[4)i 



D'F{Q,)[(f),,(f)j] = lim 



h^O h 
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(all these limits will be shown to exist). 

Moreover the diagonal terms of the covariance Cij can be written as 



where 



I {l-y{s))CovQ,{ay{s),<Pi)dM{s) 







;i.6) 



1.7) 



Jo 

is the compensated {Q^^,^ 3^) -martingale associated with the jump process o/?/[0,T]. 

2 Applications to the portfolio analysis 
2.1 Computation of large portfolio losses 

We are now going to state a definition of portfolio losses. Wlien speaking of portfolio 
losses, we mean the losses that a financial institution may suffer in a credit portfolio 
due to the default events. Many specifications may be chosen to this aim. Some 
general rules are now stated. A rather general modeling framework is to consider 
the total loss that a bank may suffer due to a risky portfolio at time t as a random 
variable defined by L^(t) = ^^Liit), where Li(t) is the loss, called marginal loss, due 
to the obligor i. Different specifications for the marginal losses Li{t) can be chosen 
accounting for heterogeneity, time dependence, interaction, macroeconomic factors 
and so on. A punctual treatment of this general modeling framework can be found 
in the book by Embrechts, Frey and McNeil [rO|. For a comparison with the most 
widely used industry examples of credit risk models see [H] or [5]. The same modeling 
insights are also developed in the most recent literature on risk management and large 
portfolio losses analysis, see [16], [I3] and [8j for different specifications. 
Here we assume that 

Li{t):=ipiyi[0,t],Xi,t), (2.1) 

where ■ ,t) is bounded and continuous in 2)[0,t] x supp(/i). In other words the 
marginal loss depends explicitly on the realization of Aj and on the history of yi. 
As a particular case of our general framework we obtain the most standard set up 
commonly used in the literature of credit risk: consider ip{yi[0,t], Aj,t) := e{Xi,t)yi(t) 
where e : ^ M"*" is a continuous function of Aj, and measures the exposure in case 
of default. Thus 

N 

LW(t) = ^e(A„t)y.(t). 

i=l 

We shall often speak of asymptotic loss or asymptotic portfolio. In this case we are 
referring at the N ^ oo case of infinite obligors. The large portfolio is intended to be 
a large but finite approximation of this asymptotic regime. 

As a consequence of the central limit theorem for the empirical measure, we obtain 
the following description for L'^it), the aggregate losses computed at time t E [0,T]: 
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Corollary 2 . Assume the reciprocity condition (R) is satisfied, and consider Li{t) 
as given in Ii2.1\) . In what follows, fory[0,t] G ©[Ojt] and A G x M, we write L{t) 
for ip{y[0,t], X,t) , and l(t) = E'^''[L(t)]. As N ^ oo the sequence 



N 



m 



converges weakly to a centered Gaussian random variable with variance V{t) where 



V{t) = E^* 



t N 2- 

m-m-P I il-yis))CovQ,{c^y{s),Lit)) dM{s, 





(2.2) 



and where M{t) has been defined in ([i. ?[) 



Proof. We apply Theorem [3] with n = 1 and 0i = = ip{y[0,t], \,t) = L(t). In this 



case 



I Hp 



N 



N 



and / 



= l{t). Notice that we can consider, 



without loss of generality, t as the final horizon of the time period, i.e. t = T. 



Remark 1 . By applying Theorem [3] with {(pj)'^^^, where (pj := (p{y[0,tj], X,tj), one 
could show that the finite dimensional distributions of the process 



N 



t>0 



converge to the ones of a Gaussian process, whose covariance can in principle be 
computed. Our Theorem [3] does not allow, however, to prove weak convergence of the 
process. 

The asymptotic expected value l(t) corresponds to the fraction of loss at time t in 
a benchmark portfolio of infinite firms. Concerning equation f l2.2p . notice that in the 
case of no interaction, (i.e., (3 = 0) we have 

Vit) = EQ'' [iLit)-lit)f]=VarQ,iLit)). 

In the case of /3 > there is a suppletive noise given by the interaction. It depends 
on the past history of the process, hence a sort of ^^memory" of the variance V{t). 



The variance given in (12. 2p involves the integral with respect to a martingale. A 
simpler form for the variance, more suitable for numerical computations, can be found 
in the following special but significant case. 

Proposition 1 . Suppose that Li(T) = e{Xi,T)yiiT) . Then as N ^ oo we have that 

'L^{T) 



N 
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where 1{T) = E'^*[e{X,T)y{T)], converges to a centered Gaussian random variable 
with variance 



V{T) = VarQXLiT))+ / {E<^'[ay{s){e{\T) - /(T))])^E«* [/?2(i-^(,))e-^+/^™Q*(^)] ds, 

Jo 

(2.3) 

where rnqXt) '■= J a'ri(t)Q^:{dri[0,T],dX'). 

Proof. We only need to show that the variance can be written in the form given in 
(12. 3p . From ( 12. 2p it is easy to see that V{T) can be written as 



V{T) = E^' 



(^L{T) - 1{T) - - y{s))E^'[ay{s){e{\T)y{T) - 1{T))] dM{t 



We now look at the expectation in the integral 

E«*M.)(e(A,TMT) - /(T))] = E«*Ms)(e(A,T) - /(T))]. 



(2.4) 



where we have used the fact that E^[y{s)y{Ty\ = E^[y{s)] for s < T and for all 
P G Ml. Substituting in ([23D we have 



V{T) = E^* 



L{T)-l{T)- / (3{l-y{s))E^'[ay{s){e{X,T)-l{T))]dM{s] 
Jo 

= E^* [{L{T) - l{T)f] - 



2(L(T)-Z(T)) / [3{l-y{s))E'i'[ay{s){e{X,T)-l{T))]dM{, 



+ E^* 



+ 



(3{l-y{s))E<^*[ay{s){e{X,T) - 1{T))] dM{, 



(2.4) 



The first expectation is the variance of L{T) computed under Q*. We show now that 
the second expectation is null. Indeed it is equal to 



T 



L{T) / p{l-y{s))E'^*[ay{s){e{X,T)-l{T))]dM{t 



- 21{T)E^' 



T 



/3(1 - y{s))EQ^[ay{s){e{X,T) - 1{T))] dM{s) 



The second term is zero since M is a (Q*; 3^)— niartingale and the argument of the 
integral is 3"^ measurable. Concerning the first one we see that 



L{T) / {l-y{s)){.)dM{s) 



e(A,T)l{^<r} 



TAt 



e(A,T) 



dM{s) 
■) dM{s 
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Where the last equahty is due to the fact that e(A,T) is 3^o measurable. Concerning 
the last term in (12.41) we now show that 



P{1 - y{s))E<^*[ay{sMX,T) - /(T))] dM{s] 



(2.5) 



[(3{1 - y{s))E<^' [ay{s){e{\ T) - 1{T))] ]\l - y{s))e-^+P^<^*^^^ds 



Indeed, being M{t) = l{r<t} — /o(l — yis))e '^~^'^"^'^*^'^^ds, it can be shown that the 
quadratic variation (M) of M, is {M)t = /^(l - y(s))e-^+'^™«*(")c/s, for all t e [0,T]. 
Thus for each progressively measurable process X = (X(t))tg[o,T] G L|^([0,T]), where 

L^([0,T]) := {X : \X{s)MM)., < oo}, 



1X1 



LU[0,T]) 



E«* r \ X{s)\^d{M)s = E'^* r \X{s)Wl~y{s))e-^^^'^Q*^'Ms. 
Jo Jo 



As a consequence, applying this last result with X = (X(t))(g[o,T] defined as 

Xit) = Pil - y(t))E«iay(t)(e(A,T) - /(T))]cit, 
by the isometry between L'^{Q,3^t,Q*) and L^^([0,T]), we have 



X{s)dM{s] 



\X\ 



L'2{Q,9t,Q*) 



Ll,{[0,T]) 



which is exactly (12.51) . Finally 



[/?(l-l/(s))E«*Ms)(e(A,T)-/(T))]]'(l-y(.))e-^+^"^'3-Wtis 



{E'^'[ay{s){e{X,T) - /(T))])' ■ E<^*[P\l - y(s))e-^+^"^'3-(^)] ds, 



and this proves (12.31) . 



In the next section we show an application of this law of large number (Theorem [2]) 
and central limit theorem (Proposition [1]). Indeed, we show the evolution of the default 
probability of certain groups of obligors and infer the corresponding probability of 
suffering large losses in a credit portfoho. 



2.2 An example with simulation results 

Consider the simplified case in which Li{T) = yi{T) and zero otherwise: the exposure 
at default is equal to one for each obligor. Moreover assume that /i = Pi5ai + (1— Pi)^A2- 
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This means that the law of the random environment A = 7) puts mass on two 

possible outcomes, that is A G {Ai; A2}. In this case the index i = 1,2 identifies a 
group of obligors with the same marginal characteristics. In other words, we split the 
portfolio in two types of obligors with different specifications. More complex specifica- 
tions could also be chosen since, as already said in the introduction, we are not forced 
to split the sample in groups with homogeneous characteristics. For sake of simplicity 
we start we an illustrative example where interesting features on the dynamic of the 
state variables can be captured. 

Recall that ai specifies the relative weight of firm i in building the aggregate variable 
rriM = aiUi (see equation (11. ip ). [3i is the parameter that measures how sensitive 
obligor i is with respect to the aggregate variable m^r, put differently it is a measure of 
the contagion effect. 7^ is the idiosyncratic term in the marginal default probability. 
As an example, we take a portfolio of = 125 obligors. This is a typical size for 
CDO's portfolios. We suppose that the portfolio consists of obligors of two types, in 
particular Ai = (4,4,3), A2 = (0.1,0.1,3). Notice that = Pi hence the reciprocity 
condition (R) applies and we are allowed to rely on Theorem [3] and Corollary [2l 
Obligors of type 1 are more sensitive to the aggregate variable, that is, their marginal 
default probability depends strongly on the default indicator of the other firms. Oblig- 
ors of type 2 are less infiuenced by the aggregate variable rriN. The idiosyncratic term 
7 is the same for each obligor. 

With this choice of the parameters, we want to stress the fact that even though the 
marginal default probabilities of the two types are rather similar for the very short 
horizon (where the impact of 7 is higher), the contagion effect becomes preeminent as 
time goes by, at least under certain specifications in the construction of the portfolio. 
To illustrate this situation, in Figure [T] (on the left), we show the dynamics of the 
marginal default probability of the two groups in two different scenarios. Scenario A 
mimics a portfolio where only 20% of the obligors are in group 1, this means that the 
proportion of firms exposed to contagion risk is lower. In scenario B the proportion is 
increased to 40%. Notice that in the second scenario the probability of default of the 
firms in the first group increases dramatically after the second year. Firms of type 2 
are less infiuenced. 

Relying on Proposition [1], we compute the corresponding excess probabilities, that 
is, the probability of suffering a loss bigger than x as a function of time in the two 
scenarios. In Figured] (on the right) we see this probabilities for x = 0.15. Notice 
that in this simple example, L^{t) counts the number of defaults up to time t, so that 
P{L^ {T)/N > 0.15) represents the probability of having at least 15% of defaults in 
the whole portfolio. Looking at the graphs, it is easy to see that in the first scenario 
the probability of having such a loss is smaller. At time t = 2.5 the ^'^^ > x) ~ 
whereas in the second scenario ^'^'^ > x) ~ 0.55. 

This simple example suggests that the contagion effect can be very significant when 
looking at the probability of suffering a certain loss in a large portfolio. This is crucial 
for risk measurements purposes and for pricing trances of CDO's. 
It may be argued that an approximation via an infinite portfolio may not reproduce 
the real situation. In Figure [2] we show a comparison between the evolution in time 
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Figure 1: Conditional probabilities of default and excess loss probabilities in a portfolio of 
N = 125 obligors where L^{t) is as defined in Section [2. 2i Scenario A represents a portfolio 
where only the 20% of obligors are of type 1 whereas in Scenario B 40% of the obligors are of 
type 1. In both Scenario A and B we have Ai = [4,4,3], A2 = [0.1,0.1,3], where Xi describes 
the idiosyncratic characteristics of obligors of type i = 1,2. 



of the default probabilities in the two groups of obligors computed under Q^, (in the 
upper part) and simulated via the real Markov process with = 125 (in the lower 
part). Here the parameters are Ai = (3,3,3), A2 = (0.1,0.1,1). These graphs show 
that the asymptotic equation is a good approximation for the Markov process even 
with A^ = 125. 

Finally we show in Figure [3] the comparison between the evolution of the aggregate 
loss (in black) evaluated as before in two different ways. In the upper part we see the 
aggregate loss computed relying on Proposition [H Below we see the plot of a trajectory 
of the Markov Process (with A^ = 125). In the same graphs we have also plotted the 
dynamics of the probability of suffering a loss over certain thresholds Ck- In other words 
we plot t ^ P{L^{t) > Ck) for = 1, 2, 3 (in this case ci = 5%, C2 = 15%, C3 = 25%). 
Those probabilities are building blocks for computing the price of trances of CDO's 
contracts. For a description of such a credit derivative see for instance [TO] . 



3 Conclusions 



We have proposed a model for credit contagion where heterogeneity and direct con- 
tagion among the firms are taken into account. We have then quantified the impact 
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Asymptotic regime: N=c« 
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Figure 2: Comparison between the dynamics of the default probabihties of two groups of 
firms computed under two different models. In the upper part we see the plot under the 
asymptotic model with infinite firms (under Q^,). In the lowest one we have implemented a 
simulation of the Markov process directly (here N = 125). The parameters are Ai = [3, 3, 3], 
A2 = [0.1,0.1,1], where Aj describes the idiosyncratic characteristics of obligors of type 
i = 1, 2. The distribution of A is in this case /i = + \6\^. 



of contagion on the losses suffered by a financial institution holding a large portfolio 
with positions issued by the firms. 

Compared to the existing literature on credit contagion, we have proposed a dy- 
namic model where it is possible to describe the evolution of the indicators of financial 
distress. In this way we are able to compute the distribution of the losses in a large 
portfolio for any time horizon T, via a suitable version of the central limit theorem. 

The peculiarity of our model is the fact that the homogeneity assumption is broken 
by introducing a random environment that makes it possible to take into account the 
idiosyncratic characteristics of the firms. One drawback of the intensity based models 
commonly proposed in the literature is the difficulty in managing large heterogeneous 
portfolios because of the presence of many obligors with different specifications. In 
this case it is common practice to assume homogeneity assumptions in order to reduce 
the complexity of the problem. A typical approach is to divide the portfolio into 
groups where the obligors may be considered exchangeable. We have shown that our 
model goes behind the identification of groups of firms that can be considered basically 
exchangeable. Despite this heterogeneity assumption our model has the advantage of 
being totally tractable: it is possible to compute in closed form the mean and the 
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Figure 3: Evolution in time of the aggregate losses (in black) computed under two different 
models: under the asymptotic model (upper) and the Markov process (lower). The param- 
eters are Ai = [2,6,4], A2 = [1,3,5], where Aj describes the idiosyncratic characteristics of 
obligors of type i = 1,2. The distribution of A is in this case fi = 0.4(5ai + O.Gdx^. 



variance of a central limit type approximation for the losses due to a large portfolios 
in a dynamic fashion. 

As an example of the general theory we have computed the default probabilities 
and different risk measures in a simple situation with only two groups of obligors. 
Moreover we have compared the numerical results obtained relying on the asymptotic 
model and on the central limit theorem (Corollary [2]) with the results obtained in a 
simulation of the underlying Markov process with finite = 125. These results show 
the goodness of the approximation and are encouraging a more involved analysis. This 
issue is left to future research: it is in fact out of the scope of this work to pursue a 
punctual calibration of the model to real data. 

A Proofs of main results 

A.l Proof of Theorem [1] 

We need to prove some technical lemmas. 

Lemma 1 . Let X be a Polish space. Let {Pn)n satisfy the LDP with rate N and good 
rate function H . Let F : X — M 6e measurable, bounded from above and continuous 
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on the set Xh '■= {x : H{x) < 00} . Then the sequence of probability measures {P^)n 
defined by 

dPL. exp(iVF(-)) 

dP^^' J^exp{NF{y))PM{dy) ^ ' ^ 

satisfies the LDP with the good rate function 



I{x) = H{x) - F{x) - m{[H{y) - F{y)]. 



In particular 



lim — loe 



exp{NF{y))P^{dy) 



X 



- mnH{y)-Fiy)] 



(A.2) 



(A.2) 



For a proof see [T9]. This is a relaxed version of the usual Varadhan's Lemma for 
tilted large deviations principles (see [20j). The statement is relaxed in the sense that 
we assume that a suitable function F : X — *• M, instead of being continuous on all its 
domain, is continuous only on a subset Xh ^ X in the following sense: 
For any sequence G X such that Xn x, where x G Xh we have F{xn) F{x). 

We point out that this is a stronger assumption than assuming continuity of the 
restriction of F on the subset X^- 



Lemma 2 . For given A G 



X 



X 



dPi 



N 



where F{Q) has been defined in U.S\} . 



iy[0,T]) = exp{NF{p^{y[0,T],X))}- 



(A.3) 



Proof. It basically follows from the Girsanov formula for point processes (See [?]). 

N „y 



dP- 



exp 



i=l 



/ {I - yi{t)) - {I - yme 



■7,+ft/a'r,{t)pjv(rf»7[0,T],dA') 



N 



1=1 



li+{f3i a'r]{t-)pN{dv[0,T],dX') 



t=Ti 



where Tj has been defined in (11.21) . 

The term in the { }-brackets is exactly -F(Q)|Q=pjv multiplied by A^. 
We now define for each Q G M and t G [0, T] 



mgit) := / a7]{t)Q{dr][0,T],dX). 



We shall often use this notation in the rest of this appendix. 



dt+ 



(A.4) 
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Lemma 3 . F(Q) is hounded on Mi and continuous on the subset JAa '■= {Q G Mi 



Proof. We rewrite F{Q) as given in fll.2p using the notation introduced in flA.4|) . 



F{Q) 



Q(dy[O.T], d\) { {I- y{t)) (1 - e-^e^'"'3W) dt 



+ y{T) 



t=Tiy[0,T]) 



(A.4) 



The argument in the { }-brackets is bounded, thus we are allowed to interchange the 
expectation with respect to Q and the time integral. 



F{Q)= I E^[l-y{t)] dt - I E<3 [e-V'"«(*) (1 -?/(t))] dt- 
Jo Jo 



- E^ [iy{T)] + E^ 



(3y{T)mQ{t- 



t=T{y[0,T]) 



We show now boundedness and continuity of F. The boundedness is easily proved 
since ?/ G {0; 1}, and the distribution of A under Q G Mj has bounded support. 
In order to prove the continuity on Ma, we consider a sequence of probabilities 
{Qn)n>o G M5 converging weakly to Q G M^- We split the proof in different steps. 
We show first that 

limE«"[/(A)y(t)] = EQ\f{\)y{t)] (A.5) 

n 

for all t G [0, T] and for any continuous function / : x R+ x M — >• M bounded on 
the support of /i. This statement is not trivial, since the projection |/[0,T] y{t) is 
not continuous in !D[0,T]. However, define for any e > the functions 



9t (e; /) := 



f{\)y{s)ds, 9tie;f):- 



t-e 



t+e 



f{\)y{s)ds; 



where we suppose that the trajectory y[0, T] can be extended to the larger interval 
[0 — e, T + e] by continuity. 

These functions are continuous in !D[0,T], bounded by ||/||oo for any e and such that 
g^{e;f) < f{\)y(t) < gf{e]f) a.s. for any t. Thus, by the Lebesgue convergence 
theorem, 

limsupi?«"[/(A)i/(t)] < limE«"[(?+(e;/)] = E%+{e-f)l Ve > 0. 

n " 

Letting e — and noticing that \mif^^Qgf{e] f) = f{\)y{t) we get 

limsup£;«"[/(A)|/(t)] < E«[/(A)y(t)]. 
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The same argument holds for (e; /); here hm£_^o 9t f) — fWuit )■ Thus 



E'^[f{X)y{t-)] < hminf E«"[/(A)y(t)] < limsnp E^-[fiX)yit)] < E«[/(A)y(t)]. 



Notice that f{X)y{t) and f{X)y{t~) may differ only on the event {y{t~) ^ y{t)}- But 
this event has measure zero for any Q e M^, since {W <S> r]){{y{t^) ^ uit)}) — 0. 
This implies that the corresponding expected values must coincide; as a consequence 
E'^[f{X)y{t)] - E^[f{X)y{t-)] = 0. We have thus proved that 



hm A)y(i)] = E'^[f{X)y{t)] for all t. 



(A.6) 



Notice that in saying that {W (E) i]) {{y (t^) ^ y{t)}) = we have used the fact that the 
distribution function of r under W ^ 7] is exponential. In particular, it is absolutely 
continuous. A similar argument shows that if Q ^ Ma then E'^"'[y{t)] converges 
pointwise in t to E'^[y(t)], for all those t such that Q{t = t) = 0. 
Taking /(A) = 1 we simply have that for all t, E'^[y{t)] is a continuous mapping in Q 
on Ma- Choosing instead /(A) = e^'^ , /(A) = —7 and /(A) = a, we prove continuity 
for E'^ [y{t) e^'^], E'^ [—'jy{T)] and mqit) respectively. 
The next step is to show that Qn ^ Q imphes that 



eQ 



/ (1 - e-V"^«(*)dt - / (l-y(i)) e-V"^«"(*)di 



converges to zero. 

We add and subtract E^^" 



EQ 



[ {l-y{t)) e-^e^"*«^*)(it - E^- [ {1 - y{t)) e-^e'^""^^'^dt 
Jo J L-'o 



+ 



where 



r ((1 - y{t)) e-^) (e^'^QW - e^"^«»(*)) dt 
Jo 



eQu 



((1 - y{t)) e"^) (e^'"^^ - e^™«"(*)) dt 



|a„| goes to zero by weak convergence. 
Concerning 6„ we see that 

rT 



■7 . \p/3mQ{t) 



\bn\< f EQ-[{l-y(t))e- 
Jo 

< I E^" [(1 - y{t)) • |e^"*«(*) - e^"'^"^*)!] dt. 
Jo 



(A.7) 
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We now show that 



We can rewrite it as 



{l-y{t)) e-^e^'™«(*) • 



^/3[mQ(t)-mQ„(t)] 



< 



j[mQ{t)-mQ„(t)] 



(*) 



for a suitable K G M"*", where we have used the fact that {1 — yit)) e '''e^"^<3(*) is 
uniformly bounded. We now look at the term into the expectation 



,f3[mQ{t)-mQ„{t)] 



1 



< K2\mQ{t) - mgjt)! 



again by uniformly boundedness, for a suitable K2. Thus 



and this converges to zero thanks to what we have shown in flA.51) . 

As a consequence, goes to zero as well, since we are allowed to interchange the 

limit and the time integral, by dominated convergence. 



It remains to show the continuity of the term 

E'^[(3y{T)m,{T-)]. (A.8) 

Indeed, take a sequence Qn ^ Q, Q ^ 3Vt^, then 

\EQ-[Py{T)mQ^ir-)]-EQ[PyiT)mQiT-)]\ 
< \E'^-[Py{T) {mQ„(r-) -mQ(r-)}]| + |i5;«"[/3y(T)mQ(r-)] -i5;«[/3|/(r)mQ(r-)]| . 

The second term goes to zero by weak convergence since the function mg is continuous. 
Concerning the first term, it is enough to show that {mQ„{t) — mQ{t)} converges to 
zero uniformly on [0,T]. To show this, we fix t G [0, T] then the following facts hold 
true: 

(a) For 61 > there exists e > such that |s — t| < e ^ Iuiq^s) — niQit)] < 61 

(b) There exists n such that \fn > n 

l^Qnit + e) - mqit + e)| < ^2 ; Imgjt - e) - nigit - e)| < ^3 

Point (a) is due to the continuity of mq^t) for Q G whereas {h) follows by the fact 
that mQ^{t) —y mQ{t) pointwise in t as shown in flA.Sp . Notice that when t + e > T or 
t — e < the inequalities in (fe) are modified appropriately without loss of generality. 
We now claim that fixing s G := [t — e, t + e] we have 

\mQ„{s) -mQ{s)\ < 61 + 62 + 65 = 6. (A.9) 
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Fix s and n and suppose mQ,^(s) < mQ^s) (the other case is treated in the same way): 

mQ{s) - mQ^{s) < mQ{t + e) - mQ{t - e) + mQ{t - e) - mQ^{t - e) < Si + S3; 

where we have used the fact that mqlt) is increasing in t for all Q G Mi. Thus 
we can extract a finite covering {Ot,,} of [0,T] where (1A.9I) holds true hence uniform 
convergence is proved. I 

Proof of TheoremUl 

We denote by the distribution of pj\j under P^, i.e. 7^ '■= P'^ ° Pn^- We now 
state a LDP for the sequence {'J'n}n- Thanks to Lemma [21 we have identified the 
Radon Nikodym derivative that relates W^^ and P^ (where PF®^ plays the role of 
the reference measure). A natural way to develop a large deviation principle is now 
to rely on Lemma [H 

Since {yi[0,T]; Xi) are i.i.d. random variables under [W ® /u)®-^, we can apply 
Sanov's Theorem (see [20]) to the sequence of measures {'Wn)n, where Wtv represents 
the law of the empirical measure in the case of independence (i.e. under {W ® n)^^). 
Hence (WAr)7v obeys a large deviation principle with rate function H{Q\W ® /i). 
Being F{Q) bounded in the weak topology and continuous on D Mh where 
3Vt^ = {Q e Mi\Q <^ W n} and = {Q E Mi : H{Q\W (g) /x) < 00}, we can rely 
on Lemma [U to conclude that the sequence (T'ly)^ obeys a large deviation principle 
with good rate function 

I{Q) = H{Q\W ® /i) - F{Q) - inf [H{R\W ® /x) - F{R)]. 



We finish the proof by showing that 

^ 0. (A.IO) 



inf [H(R\W fi) - F(R)]= lim ^ log 



e^^(«)W^(dQ) 

Ml 



The first equality is simply a consequence of Equation (IA.2p . 
We are thus left to prove that J^^^ e^^^^'>WN{dQ) = 1. Indeed, 



'Pn{-) = J /i^^(rfA)P^(p^ e ■ ) 



dP^ 



Being 'y^iJAi) = 1, the thesis follows. 
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A.2 Proof of Theorem H 

We need to define a new process and a technical lemma related to it. 

We associate with any Q e Mi the law of a time inhomogeneous Markov process on 

{0; 1} which evolves according to the following rules: 

y = Q-,y=l with intensity (^-i^^P ! a'v{t)Q(dri[o,T],d\) 
y = 1 —>■ y = with intensity 

and with |/j(0) = 1 for all i = 1, ...,N. 

We denote by P'^''^ the law of this process and by P'^ = P'^''^ (g) r/. In other words, 
pX,Q jg ^YiQ law of the Markov process on {0; 1} with initial distribution Sq and time- 
dependent generator L^'^ defined as 

Lf'^f{s) = (1 - s)e-^e'^^"'''W^('^'^l°'^]''^^)(/(l - s) - /(«))• (A.ll) 
We show now an important property of P*^. 

Lemma 4 . For every Q G Mi, we have 

I{Q) = H{Q\PQ). 

Proof. We distinguish two cases: 

Case 1 . Q : H{Q\W ® rj) <oo. We have 



By Girsanov's formula for continuous time Markov chains, we obtain 
log^^ = ^^(1 - y{t)) (l - e-^e^!Qmm4X'Wr,it)^^ dt 



l+{f3 a'r]{t-)Q{d7][0,T],dX'] 



t=T{y[0,T]) 



(A.ll) 



+ yiT) 

hence, by definition of F given in (II. 2p we have 
so that 

where the last equality follows from 

dQ dW dQ d{W0r]) dQ 



d{W®ri)dP^'Q d{W®ri) dPQ dPQ' 
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Being J log ^dQ = H{Q\P^), the thesis follows. 

Case 2 . Q : H{Q\W ®ri) = +00. In this case I{Q) = +00. 
Thus we have to check that H{Q\P^) = +00 as well. Being 

the thesis follows since J log-^^^^^^dQ = +00 being H{Q\W ^ rj) = +00 and since 
/ log j^dQ = -F{Q) which is bounded . I 

Proof of Theorem\^ 

By properness of the relative entropy {H{^\v) = ^ fi = u), from Lemma lA.llI 
we have that the equation I{Q) = is equivalent to Q = P^ ■ Suppose Q is a 
solution of this last equation. In this case Q = P^ = P^'^ ® fi, where P^'^ is the 
law of the Markov process on {0; 1} with initial distribution and time-dependent 
generator L^'*^ as defined in flA.lip . The marginals of a Markov process are solutions 
of the corresponding forward equation. This leads to the fact that qt := IlfP'^''^, the 
t— projection of P'^'*^, is a solution of qt = ^l'^ where L^'^ is the adjoint of Lf'*^: 

More specifically, when x = we have 

and when x = 1 

We now prove that qt = L'^'^q admits at most one solution for each initial condition. 
To see this, define qt{\) := P^'^iy{t) = 1). Then = G(gt(A)), where G{q) = 

g-7g/3/a'g(A')M('iV)(i _ g(^x)). Notice that qt{X) e L\fi) and G(-) is a locally Lipschitz 
operator on a Banach space. Thus = ^(^^(A)) has at most one solution in [0, T], 
for a given initial condition (see [2|, Theorem VII. 3). 

Since P'^ is totally determined by the flow qt, it follows that equation Q = P'^ has at 
most one solution. The existence of a solution follows from the fact that I{Q) is the 
rate function of a LDP, and therefore must have at least one zero. By what shown in 
f lA.13p . gt(A) solves fll.3l) . Hence turns out to be the unique solution of the fixed 
point argument Q = P^ . Moreover, it satisfies all the conditions of Theorem [21 I 

A.3 Proof of Theorem [3] 

The key technical tool for the proof of Theorem [3] is the following result due to 
Bolthausen p]. 
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Theorem 4 . Let {B, \\-\\) be a real separable Banach space. Let {Zk)k>i be a sequence 
of B -valued, i.i.d. random variables, defined on the probability space {Q,A,F) and 
denote by w their common law. Define Xn := ■^k o,nd consider a continuous 

map \1/ : i? — > R. Suppose that the following conditions are satisfied: 



(B.l) J exp{r\x\)w{dx) < oo for all r G M. 

(B.2) For any x E B, "ifi^x) < Ci + C2||a;||, for some Ci, C2 > 0. Moreover, \1/ is three 
times continuously Frechet differentiable. 

(B.3) Define, for h E B' (the topological dual of B), A{h) := J e'^^^^w{dz), and for 
x & B, K*[x) := supf^f^^,[h{x) — A{h)]. Assume that there exists a unique z* E B 
such that A*{z*) - ^{z*) = inf^eB[A*(z) - ^{z)]. 

(B.4) Define the probability p on B by ^ = ^'"'^^^"^ for a suitable normalizing factor 
c. This probability is well defined and J zp{dz) = z* . Let p^, denote the centered 
version of p, i.e., p^ = p o 6~} , where 6a '■ B ^ B is defined by 9a{x) = x — a. 
For h E B' define h E B by h = J zh{z)p^:{dz) . Then we assume that for every 
h e B' such that h^O 

j h^{z)p,{dz) - D^^{z*)[h, h] > 0. 

(B.5) B is a Banach space of type 2. 
Now, letting TTjy be the probability on {fl,A) given by 



(A.14) 



where S is linear, continuous and bounded on the support of the law of X^, uniformly 
in N . Then, for every hi, hn G B' , the Tt^-law of the n— dimensional vector 



N{K{XN)-h,{z,)YU 

converges weakly, as N ^ 00, to the law of a centered Gaussian vector with covariance 
matrix C G M"^", such that for i,j = 1, ...,n 



{e)ij= hi{z)h,{z)p,{dz)-D'm{z,)[hi,h,]. (A.15) 



Remark 2 



i. The Theorem in [3j is stated for S = 0. The same proof apphes with our assump- 
tions on S without changes. It is hkely that these assumptions can be weakened 
considerably. 



1 A Banach space Bis said to be of type 2 if ^ (jl^^y Here^^j-^) ^ {(^.^j ^ ^00 . ^J|a;j2 < 

00} and C{B) = {(x„) G : Ylij ^j^j converges in probability} where (e„) is a BernouUi sequence, 
i.e., a sequence of independent random variables such that P(e„ — ±1) — ^. For more details see [5]. 



23 



ii. The Theorem in [3] contains a stronger statement that the one given here. Indeed, 
it is shown that the field (/^(Xtv) — h{z^))^^Qi converges weakly to a Gaussian 
field, while we only stated convergence for finite dimensional distributions. This 
is all we need to prove Theorem [3J Our proof has not allowed to take advantage 
of the full strength of the Theorem in [3]. 

The "natural" space for the Central Limit Theorem in Theorem |3] is the set M of 
signed measures, which is not a Banach space. To apply Theorem IH we need to map 
M to a Banach space of type 2. 

Lemma 5 . The following properties hold true under the reciprocity condition (R): 

i) There exists a Banach space of type 2 {B, || ■ ||), a linear map T : M. ^ B, 
continuous on the set {Q : Q{t = T) =0}. Moreover there exist two continuous 
maps \E', S : i? — >■ M, where \E' is bounded and three times Frechet differentiable 
and E is linear, such that 



dP, 



dW®^ 



A 

^ -exp<'A^ 



"^{TiPN)) + 



N 



a.s. (A.16) 



ii) For any vector $ = e Cfe"" there exist h = {hi, hn) C B' such that 

{hi o T){Q) = J ^idQ, where B' stands for the topological dual of B. 

Proof. The first step consists in giving an alternative expression for F{Q) given in 
(11. 2p . Look at the term 



t=T{y[0,T]) 



Q{dy[Q,T],dX) \y{T) (^/? j Q{dr][0,T],dX')a'r]{t- 
using the reciprocity condition (R), we obtain 

b [ Q{dy[0,T],dX) L{T)a [ Q{dr][0,T],dX')a'r]{t- 



t=T{y[0,T])^ 

= b j Q(rf?/[0,T],rfA)(5(c/r7[0,r],c/A')aa'l{r(j/[o,T])<T}l{T(r,[o,T])<r(?;[o,r])} = (*) 
Note that l{r(,?[o,T])<r(y[o,T])} = l{r(,,[o,T])<T} uuless r(?/[0,T]) < r(r7[0,T]) < T. Thus 

{^) = b j Q(c?y[0,T],dA)Q(ci77[0,r],c/A')aa'l{r(y[0,T])<T}l{r(,7[0,T])<T} 

-b j (5(dy[0,T],ciA)Q(rfr/[0,T],rfA')aa'l{r(,?[o,T])<T}l{r(y[o,T])<r(,?[o,T])} 

= b j Q((i?/[0,r],(iA)(5(rf?7[0,T],rfA')aa'l{rfe[o,T])<r}l{r(,,[o,T])<r} 
-b j Q{dy[Q,T],d\)Q{dri[Q,T],d>^)aa'l{r(^r,[o,T])<T}'^{T{y{o,T]^^^ 

-b Q{dy[0,T],dX)Q{dr][0,T],dX')aa'l{r{r,[o,T])<T}^{T{y[o,T])=T{v[o,T])}- 
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Therefore 



b / Q{dy[0,T],dX)L(T)a / Q{dr][0,T], dX')a'r]{t\ ^_ 

= 2 / '5(c?y[0,T],dA)(5((i?7[0,T],(iA')aa'l{^(j/[o,T])<T}l{r(r,[o,r])<T} 

- 2 / <5(c?2/[0,T],rfA)Q(dr7[0,T],rfA')aa'l{r(r,[0,T])<T}l{r{j/[0,T])=r(r,[0,T])} 



(5((i?/[0, T], rfA)al{^(y[o,T])<T} 



H 2 T 

- E 



te[o,T] 



(5(rf?/[0, T], rfA)al{^(j,[o,T])=t} 



Thus, defining 



Fi(Q) := / Q{dy[0,T],dX) 



(1 _ y^t)) (l - e-^e^/«('^''[0'^]''^^')-'''W) - 7y(T) 
6 



+ 



n 2 



g(ci2/[0,T],dA)a2/(T) 



and 



HQ) -=-2^1 J Qidy['^,T],dX)al{riy[o,T])- 



tG[0,T] 



(A. 17) 



(A.18) 



we have that F{Q) = Fi{Q) + F2{Q). Lemma [2] thus holds also after replacing F by 
Fi + F2. Let M > be a constant such that, under rj, the random parameters a and 
7 have absolute value less that M/2. Now we define the following maps: 



Ti : 



M 



L2[0,T] 



Q ^ UQidy[0,T],dX)[l-y{t)] 



M 
Q 



L2([0,T] X N) 
CiM)jQ{dy[0,T],dX) 



T3: 
T4: 



T5: 



M - 


- L2[0,T] 


g H 


/Q(rfy[0,T],dA)^y(t) 


M - 




Q H 


jQ{dy[0,T],dX)^y{T) 


M - 


R 


Q H 


/g(t;y[0,T],rfA)^t/(T) 


M - 


> M 


Q ^ 


. /g(t^l/[0,r],f/A)a2^(T) 
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where C(M) is some positive constant such that C(M)e~'^e^^ ^ | ^ almost surely. 
Note that, for Q G Mi and i = 1,2,..., 5, we have that < |. Now, let 

5f : M ^ M be a C°° function such that g{x) = x for |a;| < 1/2, g{x) = for |a;| > 3/4 
and lls'lloo < 3/4. For 

Z = {Zi, Z2, Z3, Z4, Z5, Zq) G L'^[0,T] X L'^{[0,T] X N) X L'^[0,T] X R X R X R 

we set 

vi/(z):=2y^ 9i^imt--^^J2j^ 9{{z,{t,n))g{z,{t)rdt-Mg{z,) + ^g\z,) 

(A.19) 

and 

:= -^-z,. (A.20) 
We now claim that, for Q G Mi and setting T = (Ti, T2, T3, T4, T5, Te), 

Fi(g) = vi/(T(g)). (A.21) 

Moreover 

F2(P7V) = H^®/x-a.s. (A.22) 

so that (]A.16P holds. Equation (]A.2ip is straightforward; (1A.22P follows since 

t \ i / it i 

where the one to last equality follows since simultaneous jumps may happen only with 
zero (ly probability. We thus have that F2{pn) = J ce'^y{T)dpN and the 

claim follows by definition of S. 

Now set 

B := L'^[0,T] X L2([0,T] x N) x L'^[0,T] x M x M x M x M". 

Clearly i? is a Hilbert space (hence a Banach space of type 2), and the maps \E', S are 
trivially extended to B. Moreover, the map T can be completed to a 5-valued map 
by letting, for i = 1, 2, . . . , n. 



where $ = ($1, $„) G 6^" is given. 

To complete the proof of part i) of Lemma one has to show the desired regularity 
of ^ and S. The only nontrivial fact is to show regularity of the term 

^ f g{{z2{t,n))g{z,{t)Ydt. 

n 
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However, the fact that ||5'||oo < 3/4 allows to control the tails of the sum above; 
continuity and Frechet differentiability of any order is obtained by standard estimates, 
The details are omitted. 

Finally, to prove part ii), for B 3 {zi, . . . , zq, . . . , ^e+n), it is enough to define for 
i = 1,2, . . . ,n 

hi{z) = ZQ+i. 



Proof of TheoremlBi 

Having identified a suitable Banach space. Theorem [3] immediately follows from The- 
orem H] applied to the sequence Zi = T{6{y-[o^T],\}) taking values on {B, \\ ■ ||). Notice 
that in our setting n = {D[0, T] x R'^)^ and P = (V (g) /i)^. Theorem [3] is guaranteed 
by the following three facts: 

1. P/v = ttat, where njy is the probability appearing in Theorem 01 

2. (/ ^dPN - J $iC^Q*)r=i = (HXn) - /i.(^*))r=i; 

3. /(0,-0*)(0,-0*)dQ,-D2F(g,) = Jh,{z)h^{z)p,{dz)-D^^{z,)[k,h^]. 

Point 1. follows from the definition of vttv and from equations ( 1A.21[) and ( ]A.22l) . Point 
2. is a consequence of the fact that 2* = T{Q^) and hioT{Q^) = f ^idQ^. Point 3. will 
be proved in details in Lemma [H] (see in particular equation ( 1A.35I) ). An immediate 
application of equations (1A.35P and (1A.30P finally guarantees the validity of (11.61) . 
Assuming point 3., it remains to show the validity of the central limit theorem in B. 
In other words we need to check the five assumptions of Theorem HI (-B.1), {B.2) 
and (-B.5) are easy to see. {B.3) and (BA) are not straightforward. The rest of this 
section is devoted to the proof that these two assumptions are satisfied. 
We begin to prove {B.3). We define two sequences of measures on B as follows: 

Pn{-) = -Pn o T-\-) wn{-) = o T-\-). 



From flA.16p it can be shown that 



dw]\f 



(A.23) 



for ^ and E as defined in flAlOj) and ([OOj) . 



By the contraction principle (see Theorem 4.2.1 in [9J), the sequence {pn)n satisfies a 
LDP with the good rate function J{z) = miQ=T-^(z) HQ)- Being the unique zero 
for /, J has a unique zero = T{Q^). 

A LDP for the sequence {pn)n can be obtained in an alternative way. Indeed, we 
notice that wn is the law of the random variables 



where Zi are i.i.d. i?— valued random variables with law w. Thus we have that {wn)n 
satisfies a (weak) LDP with rate function A*, with A*{z) := sup^^^, {ip{z) — A{ip)} 
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and A(y9) := In J e*^*^^^ w{dz). Thus, applying Varadhan's Lemma, {pn)n satisfies a 
(weak) LDP with rate function A*{z) — "^{z). Since the rate function is unique, it 
follows that 

J{z) = A*(z) -^(z). 
Having proved already that J{z) has a unique zero, the proof of (B3) is completed. 

We are thus left to show (-B.4): for each X E B' such that h = J zh{z)p^{dz) ^ 
we have 

j h\z)p4dz) - D^<il{z*)[h, h] > 0; (A.24) 
where p and are defined in Theorem HI 

This proof is rather technical and long. We divide it into three steps. We first show 
that the measure p such that ^ = e^"^^^*^ is exactly the law of the random variable 
T(5{y[o_T],A}) induced by Q*. This argument is then used in the second step to ensure 
the positivity of a suitable functional : 6;, x C;, M. In the last part we see how to 
relate "K to assumption (-B.4). 



Step 1: The key result of this first step is given in Lemma [6] below. We look at the 
measure p on B, defined by 

{z) = e^*(^*)W , with z, = T{Q,) 



dw 

where, as already seen, w represents the law of T{6{y[o^T],x}) induced by ® /i. 
We shall prove in Lemma [6] that p is the law of T{6{y[o^T],\}) induced by Q*. 

Lemma 6 . The measure p is the law of T{6{y[o^T],x}) induced by Q^. 

Proof. We first prove the following two claims 

i) 

DFmiS^y^onx}] = log ^^^^^^^ (l/[0,T],A), (A.24) 
for W (g) /i-almost all {y[0, T], A). 

ii) 

DF2{Q)[r]=0 (A.25) 
for all Q G Ml such that J a'l{r{y[o,T])=t}dQ = 0, where F2 is defined in ( jA.181) . 

To prove the claim we need to compute DF{Q^)[6{y[Q^T],x}], i-e. the Frechet derivative 
of the function F at in the direction (5{j/[o,t],a}- An explicit computation reveals 
that for Q G Mi and r G M, DF{Q)[r] is well defined and in particular 

n^(n^\^ r ^(Q + hr) - FjQ) 
DFmr] = lim = 

dQ (1 - y(t))m,(t) e-^+'^'"«Wrft + [ dr [ (1 - ?/(t))(l - e-^^^"'Q^'^)dt+ 
+ 







J dr [-y{T){^ - f3mQ{T-))] + J dQ [y(T)m,(r-)] 
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where we have put as usual mp(t) = f a'ri(t)p{drj[0,T], da') for p G M, t G [0,T]. 
We now compute DF{Q^)[6^y[o^T],x}]- Notice that m^j^^^ (t) = ay(t) = for all 
t < T{y[0,T]) so that 

DFm[6{y[o,nx}] = f{l - y{tm - e-^+^'"«*W)rft - [y{T){^ - /5mQ.(r-))] . 
Jo 

(A.26) 

By virtue of Girsanov's Formula for Markov chains it can be seen that 



where P*^ is the law of the Markov process with generator given in flA.lip . flA.241) 
thus follows since P^' = as shown in the proof of Theorem [2l Concerning (u) , 
notice that 



F2{Q + hr) = -^J2^j »'Mr{ylo,T])=t}d{Q + hr} 

= -\^(^j (^"^{r{y[0,T\)=t}dQ^ ~ J] (/ «'l{T(y[0,T])=t}C?Q j (^''^{r{y[0,T])=t}dr 

When J a'l{r{y[Q,T])=i}dQ = 0, the first two terms of the last expression vanish. More- 
over under the same hypothesis F2{Q) = 0. Hence 

DF2{Q)[r] = lim - [F^iQ + hr) - F^iQ)] = E (/ ^'Mriy[o,T])=t}drj = 0. 

(A.27) 

Notice that in writing the latter equality we have implicitly used the fact that 
J2t (/ '^''^{T(y[o,T])=t}dr) < oo. This is true since for any r G M 



< X] a'l{r(y[o,T])=t}cirj < || a\ 



r\rpy < OO 



where ||a;|| stands for the supremum of a in the support of and | r\Tv denotes the 
total variation of r. 

As a corollary of claim [ii) above we see that 



h~*0 h 

lim fiW.+'"-)-fiW.) = r>F,{Q.)\r] = Z)F(Q,)|, 
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Here we have used (1A.27|) . the fact that J (y'l{T(y[o,T])=t}dQ* = since ^ ® 
and equation (lA.2ip . 

Back to the statement of the lemma, we see that for h measurable and bounded 



h{z)p{dz) = / /i(z)e^*(^*)Ww;(rf2) 



where in the last equality we have used (]A.24p . ■ 
Step 2: The key result of the second step is given in Lemma [7] below. It involves the 



measures ^ and $* defined in f ll.Sp . First of all, it is not difficult to show that $ is 
absolutely continuous w.r.t. Q* i^i particular 



dQ^ 



(A.28) 



Indeed, observe that, given $ as in (11.51) . we have 



Mo 



R{S) j MR 



i^*{dR) 



T>[0,T]xSupp{fi) 



L{(j/[0,T],A)eS} 



Q*{S)) 



$ d6. 



{y[o,T],x} 



MQ, dQ 



[(l|(,[o,T],A)e5} - Q*{S)) ■ (<f (y[0,T], A) - $*)] dQ, 

'D[0,T]xSupp(^i) 

for any S G T)[0,T] x Supp{fi). The second equality follows since z/* is the law of the 
random variable 5{y[o,T],A} — Q* induced by Q,. 

Notice that Q,{S) /x)[o T]x5«pp(/i)('^ ~ ^*)dQ, = 0, being $* the expectation under 
of $(■). Hence 

<l(5) = / (<l>-<l>*)rfg, 
Js 



and (1X28]) follows. 

Lemma 7 . Given $i and $2 in Gb, let 

$2) := CovQ^i^i, $2) - D^F{Q,)[^u l>2] 
where Cowq,(<I>i,$2) := - $1)($2 - '^DdQ,. Then 

3<{^, $) > 0, for all $ such that $ 7^ 0. 



(A.29) 
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Proof. A tedious but straightforward computations provides the second order deriva- 
tive of F: 



1 - y(t))/32m^^(t)m4 rt)e-^+^"^«*Wdt 
- ($, - $*) / (1 - y(t))/5m|,.(t)e-^+'3™«*Wdt 



'y{T)Pm^{T-)\+E^^ [y{T)Pm^^{T-)\ . 



Notice that we have written (3 instead of ha: the reciprocity condition is not necessary 
in this calculation. We now show that IK($, $) is the expected value of a square. 
Indeed 



= ^Q.[($_$*)2] + ^Q* j 

+ 2(<l> - $*) / (1 - ?/(t))/5m^(t)e-^+'^"^«*Wdt 
The latter expectation can be rewritten as 



2E^ [yiT)P^ 



[y{T)f3' 



($_$*) / {l-y(t))(3m^{t-)dN{t) 



where we have used flA.281) and where (A^(t))ig[o,T] defined by N{t) := l{T>t}, is the 
Poisson process with intensity Jg*(l — y{s))e~'^^^'^'^*'^'^''ds. 

Recall that M{t) = N{t) - /^(l - y{s))e-^+'^"'Q*^''^ds defined in (O) is nothing but 
its compensated Q*— martingale. Hence 

r-T 



$) = E^* [($ - <l>*)2] + E^* 



'l~y{t))P'[m^{t)]'e-^^^'^^'^'^dt 



2($_$*) / (l-y(t))/3m|,(t)ciM(t) 



By the isometry property of square integrable martingales (and relying on the same 
argument used to prove (12.51) ). we have 



E^' / (l-|/(t))/?2,ri^(t)V^+^"«*Wdt 
Jo 

Hence 

J{(<l', $) = E^* (($-$* 



{l-y{t))(3m^{t)dM{t) 



T X 2' 

{l-y{t))Pm^(t)dM{t) 



(A.30) 
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$) is thus the expected value of a square, hence it cannot be negative. For this 
reason, we simply need to prove that it is non-zero. Without loss of generality we take 
$* = 0. Suppose by way of contradiction that ^K($, $) = 0. Then necessarily 

$(y[0,T],A)-^ il-yit))Pm^{s)dM{s)^ =0, a.s. 
Using the fact that 

m^is) = J ay{s) ^(^^[O, T], rfA) = j ay{s) <I>(2/[0,T],A) Q,{dy[Q,T],d\), 

where the last equality follows since = $. We rewrite the expression above as 



ay{sMy[0,TlX)dQ, 



$(2/[0,T],A)= / il-y{t))(3 
Jo 

On the other hand, define $j = -E*^* [^jJ't], where 

S't = a{ys ■.0<s<t; A}. 

Notice that 

ymi-)dQ, = E«*[ay(t)$(-)] = E^'[ayit)E'^''m-)\3^t 



dM{s) , Q^-a.s. (A.31) 



Taking the conditional expectation in flA.311) . we obtain 



l-y{t))[3 (^j ay{s)^sdQ}j dM{s] 

.] 



1 - y{s))(3 ay{s)(!>sdQ,j dM{s), Q, a.s. 
We now take the L^-norm in both sides. For all t G [0, T] we have 



ay{t)^tdQ* 



l-yis))p{ / ayis)^sdQA dMis) 



j .y,s)..dQ)\-^^^-^^^^^ds 



Notice that (/ ay{s)^sdQ,f < (/ ^ "^^sdQ* < J ^jdQ, = 

||a;|p||$( ||^2(Q^'), where t > s. The first inequality follows since y G {0; 1}; the second 
one is trivial and the latter one is due to the fact that is a submartingale and 

thus its expected value is an increasing function of time. Then 



1$ 



t\\L\Q,) 



< WafE^* 



Jo 



< te ' ||$i||i2(Q.) 



32 



where < e < oo is a constant such that ||Q;|p||/3|p£"3* [e-T+^'"«*('')] < e"^ As a 
consequence, $s = 0, a.s. for s G [0,e). 

(2) 

This argument can be iterated defining $J := ^t+e- The same argument shows that 
$s = 0, a.s. for s G [0,e); hence $t = 0, Q^, a.s. for s e [0,2e). Eventually we 
extend the statement to s e [0,T]. Being $r = ^5 we would have $ = and this 
gives a contradiction. Hence the thesis follows. ■ 

Step 3: Consider Ai, A2 G B'. Since Aj o T, for i = 1,2, are in the topological dual of 
M, there exist $1, $2 e Cft such that A, o T{Q) = J ^idQ. We define 

= I ^i{z)\2{z)p^{dz) and A, = z\i{z)p^{dz); i = 1,2 



where we recall that defined in (-B.4) of Theorem HI is the centered version of the 
law of T{6{y[o^T],x}) induced by Q*. Then the following result holds true 

Lemma 8 . 

CoVp,{Xi,X2) = Cot;Q,($i,<l>2); 

ii) For Xi, i = 1,2 we have 

QoVpSKk) - D^^{z*)[~k,~Xi] > 0. 

Proof. Point (z). By the definition of and Ai, z = 1, 2 we see that 

C0Vp^{Xi, X2) = J [T6+i{S{y[o^T],X}) - ^6+l(Q*)][T'6+2(%[0,T],A}) " ^"6+2 = 

where we have used the fact that Xi o T{Q) = T^+i^Q) = f ^idQ. 
Concerning the second statement, we first prove the following claim. 

Xi=T{^i) ; 1 = 1,2. (A.32) 

To show the validity of flA.321) . we use the following two facts: 

Xi = EQ*{[T(5{,[o,T],A})-r(g,)][$i(|/[0,T],A)-$*]}; 
1), = EQ'{[6{y[o,T]M-Q*miyM,X)-^*]}. 

The former follows by definition of p*, A and T^^ilQ), whereas the latter is a conse- 
quence of the definition of $ given in fll.51) . 

(1A.32P is a consequence of the fact that T is both linear and continuous, hence we are 
allowed to interchange the T operator with the expectation. 
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Having proved (1A.32P , we compute the second order Frechet derivatives on the function 
\1/ as follows. 

c'^*(..)|A,Aj = im,^i<£i±^MM^££MM. (A.33) 

k^O k 



Notice that, by the linearity of T and by flA.32p . we have that 



z, + kX2 = T{Q, + k^2) , z, = T{Q,). 

Thus 



k^O k k^O k 

We now claim that 



k-*0 k k^O k 

(A.34) 

By ( 1X25!) we see that DF2iQ*)[-] = since < {W^t]). Moreover both DF2{^i)[-] = 
and DF2{Q^: + k^i)[-] = since $j is absolutely continuous with respect to Q*. This 
proves ( lA.34p . Finally we use the fact that F is Frechet differentiable 

k^O k 

We have thus proved that D^^{z^)[Xi, A2] = D^F{Q^)[^i, ^2]- 
Concerning point (ii), we notice that 

CoVp^iK, Xi) - D^^{z,)[Xi, Xi] = CovQ,{^i, $,) - D^F{Q,)[^i, |),] = 

(A.35) 

where J{ has been defined in Equation ( ]A.29p . Hence by Proposition [7] the positivity 
condition is ensured and the thesis follows. ■ 

By virtue of Lemma [HI for any X E B' such that A 7^ 0, ( ]A.24p holds true. As a 
consequence, assumption (-B.4) is ensured and thus Theorem [3] is proved. I 
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